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Abstract 



We consider simple examples of self-organized critical systems on one-dimensional 
superlattices without local particle conservation laws. The set of all recurrence states 
are also found in these examples using a method similar to the burning algorithm. 
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I. INTRODUCTION 



In 1987, Bak et al. proposed an interesting idea called self-organized criticality, suggesting that 
many physical systems could evolve under their own dynamics without any fine tuning parameters 
to states without any characteristic time and length scales JIJ. They illustrated the idea using a 
simple cellular automaton model. Later on, their automaton model is proved to be commutative 
in the sense that the order of "particle addition operations" does not affect the outcome of the 
final state 0. Now, this model is widely known as the Abelian Sandpile Model (ASM). In addition 
to the cellular automaton models, various kinds of continuous and lattice-continuous models have 
shown to exhibit self-organized criticality. 

In particular, the self-organized critical model of dissipative transport using stochastic partial 
differential equation || has received much attention recently. This model proves that the existence 
of a local (particle) conservation law is not a necessary condition for the exhibition of self-organized 
criticality. Further discussions on the importance of particle conservation law can be found elsewhere 
Furthermore, the existence of local particle conservation law is not an essential feature for 
cellular automaton models (i.e. models on finite grid points and each grid point allows only a finite 
number of states) as well. The well-known forest fire model is an example of this kind ||. 

In this paper, we study two new examples of Abelian Sandpile Models on one- dimensional 
superlattices — the simplest possible kind of self-organized critical model without local particle 
conservation laws. The examples we provide here can be generalized easily to higher dimensions 
although exact calculation of the scaling exponents and the recurrence phase space configurations 
is much more difficult. Studies of similar one-dimensional models can be found elsewhere ||. A 
useful method, inspired by the idea of the burning algorithm @, to find the recurrence phase space 
which constitute an important step in showing the criticality of our model is also introduced. This 
method works well for systems of low spatial dimensions. 

II. EXAMPLE 1 

Consider a collection of 2N sites labeled from 1 to 2N and we associate an integer hi called 
the local height to each site. Using the same rules for the Abelian Sandpile, sites with local height 
greater than is said to be unstable and particle redistribution occurs as follows: (a) for any odd- 
numbered site 2i+l, two particles will be lost in the next timestep. Each of its two nearest neighbors, 
i.e. sites 2i and 2i + 2, receives two particles in the next timestep. Clearly two extra particles are 
created in the particle redistribution (or toppling) process, (b) for any even-numbered site 2i, four 
particles will be lost in the next timestep. Each of its two nearest neighbors, i.e. sites 2i — 1 and 
22 + 1, receives a single particle. Thus, two particles are dissipated in the process. The odd and 
even-numbered sites are called creative and dissipative respectively. Open boundary conditions (i.e. 
particles are allowed to flow out of the system from both ends) are used. Therefore, the toppling 
matrix |2j of the system is given by 
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A = A< 2JV > 



2 -2 
-1 4 -1 
-2 2 -2 
-1 4 -1 



(1) 



where the superscript (2N) denotes the total number of sites in the system and hence the size of 
the matrix. Clearly we have constructed a one-dimensional superlattice model with equal number 
of dissipative and creative sites. 

First we find the total number of recurrence states on the system which equals det A p[ . By 
some elementary row and column transformations, it is easy to show that 



det A< 2fc > = 4 det A (2/c-2) - 4 det A^" 4 ) 
for k > 3. Since det A^ 2 ^ = 6 and det A^ 4 ^ = 20, we conclude that 

det A^ = (2N + 1) 2 N . 



(2) 



(3) 



Now we want to find all the (2N +1)2^ recurrence states of the system using an idea inspired 
by the burning algorithm p] . We represent system configurations by row vectors of length 2N. It 
has been shown that a system configuration a = (oti, . . . , cv2n) is in the set of all recurrence states 
Q if and only if we can find an unstable state (3 = (j3i, . . . , /?2jv) which topples to a with all the sites 
topple at least once in the process |y|. Clearly, $ = ctj + J2j n j^ji f° r some rij G Z + . To test 
if a particular state a is recurrence, we choose a set of G Z + with Oj = Y^j n j^ji — f° r an 
The existence of such set of n, has been proved in Proposition 1 of Reference [1(J . Nevertheless the 
choice of n, is not unique. We add particles to site % for all i at the same time when the system 
is in configuration a. Then the resultant configuration after toppling equals a if and only if a is a 
recurrence state |9j,PH. 

In the present case, we choose rii = 1 for all i, which is equivalent to adding a single particle to 
sites 1 and two particles to site 2N at the same time. Since at most An particles are removed from 
site i whenever it is unstable each time, the possible local heights an odd (even) site can be when 
the system reaches its recurrence phase space are -1 and (-3, -2, -1 and 0). An odd (even) site 
is called "absorbing" if and only if its local height equals -1 (-2 or -3). Obviously, these are sites 
which can "absorb" the particles coming from their neighbors during toppling for exactly one time. 

Claim: any system configuration with more than one absorbing site is not a recurrence configuration. 
Proof: Suppose a has more than one absorbing site, and the left most and right most absorbing 
sites are denoted by / and r respectively (i.e., site i is not absorbing if % < I or % > r). Upon addition 
of a particle to site 1 and two particles to site 2N, it is easy to verify that exactly one toppling will 
occur in sites 1,2,...,/ — 2,1 — l,r + l,r + 2,..., 2N. After that, the avalanche stops because both 
I and r "absorb" the incoming particles and prevent further toppling. Since / < r (or else a has 
only one absorbing site), the system does not return to a after the avalanche. Thus a ^ Q. □ 

The remaining possible recurrence state configurations are those a with at most one absorbing 
site. The total number of such states, T, is given by: 
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N N 



T = 2 7V + ^2 7V + ^2 2^ = (2iV + l) 2 W (4) 

2=1 1=1 

where the first term is the total number of stable configurations without any absorbing site, the 
second (the third) terms are the numbers of possible recurrence configurations with exactly one 
odd and no even (one even and no odd) absorbing site. Since T = det A, we conclude that stable 
configurations with at most one absorbing site are the only elements in the recurrence phase space Q. 
The above method of finding recurrence configurations is effective whenever the spatial dimension 
of the system is low. 

Having explicitly find out all the elements in Q, we can proceed to show that the system is 
indeed self-organized critical. We define the avalanche size s to be the total number of toppling 
occur during an avalanche (i.e., sites with multiple toppling are counted multiple number of times). 
Direct calculation tells us that if we add a single particle to site i on a recurrence system configuration 
without any absorbing site, 

i(2JV+l-i) if hi = Q . s 

otherwise 

Similarly, if the particle is added to a recurrence configuration with an absorbing site at k, 

if i = k or hi < 

(2N + 1 — — k) if i > k and h { = . (5b) 
iik —%) if i < k and hi = 

So under a uniform particle addition and in the N — > oo limit, the distribution of avalanche size 
can be well approximated by 

\ if 8 = 

D(s) = { IW^-w] ifO<,<iV 2 . (6) 
otherwise 

Thus if we Fourier transform D(s), 1/f 2 scaling in avalanche size s is observed. Therefore this 
model is self-organized critical although its scaling exponent is trivial. 

To some extend, the above example does not completely demonstrate that particle conversation 
law is not a necessary condition for the exhibition of self-organized criticality in cellular automaton 
models. If we rescale all the even-number sites by hu — >/i2i/2, then the toppling matrix becomes 
that of the one- dimensional Abelian Sandpile Model. In fact, this example is equivalent to the one- 
dimensional ASM which allows half-integral local heights for all the even-numbered sites. Besides, 
the even-numbered sites will receive only one half of a unit of particle each time when something is 
dropped onto them. The distribution of avalanche size D(x) can be calculated using the piecewise 
linear relationship found by Chau and Ho \T2\ and it turns out to be the same as those given by 
Eqs. (|) and (|). 
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III. EXAMPLE 2 



We now provide another example on one-dimensional superlattice which is not equivalent to any 
"ordinary" one-dimensional sandpile models, whose toppling rules are translational invariant at all 
sites except possibly at system boundaries, under local rescaling. More precisely, this model is not 
equivalent to any one- dimensional ASM whose toppling matrix is of Toeplitz form. Again, we apply 
the rules of Abelian Sandpile to a collection of 2N sites labeled from 1 to 2N . The toppling rules 
are given by: (a) for any odd-numbered site 2i + l, it will lose four particles in the next timestep. 
Two of them are delivered to site 2i + 3, and one of them to site 2i + A while the remaining one is 
dissipated in the process. Thus, odd-numbered sites are dissipative. (b) for any even-numbered site 
2i, it will lose two particles in the next timestep. Two of them are transported to site 2i — 3, one of 
them to site 2i — 2. Thus, even-numbered sites are creative. The toppling matrix of the system is 
given by 



A EE A< 2JV > = 



4 





-2 


-1 





2 










4 


-2 -1 


-2 


-1 





2 



(7) 



Using similar idea as in Example 1, it is easy to show that 

det A< 2JV > = (N + 1) 4". 



(8) 



We use the same technique before to find all the (N + 1)4^ recurrence states of the system. In this 
case, we choose rii = 1 for all i, which is equivalent to adding four particles to sites 1 and 2N — 1, 
together with two particles to sites 2 and 2N all at the same time. Since at most An particles are 
removed from site i whenever it becomes unstable each time, the possible states an even (odd) site 
can be when the system reaches the recurrence phase space are -1 and (-3, -2, -1 and 0). For any 
stable system configuration a = (hi), we define 



u, 



i if hi, hs, . . . , hi 
N+l if h ± ,h 3 , . . . , h 2N -i 



3 > -1 and h 2 i-i < -1 
> -1 



and 



lr 



% if h 2i+2 , h 2 i+4, • • • , h 2 N 
if h 2 , h A ,..., h 2N = 



and h 2 i = —1 



(9a) 



(9b) 



For the odd-numbered sites, exactly two particles are received per toppling in the particles flows 
into them; while for the even-numbered sites, the only one particle is received per toppling. So 
odd-numbered sites with local height equals -3 or -2 (or even-numbered sites with local height 
equals -1) can "absorb" the particles coming from their neighbors during toppling for exactly once; 
while odd-numbered sites with local height equals -1 or (or even-numbered sites with local height 
equals 0) becomes unstable whenever particles flows into them during an avalanche. Thus 2u m i n —\ 
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is the minimum odd-numbered site (while 2l max is the maximum even-numbered site) in our finite 
superlattice which can "absorb the stress" whenever particles topple onto it for exactly one time 
during an avalanche. 

Claim: any system configuration with u m i n < l max is not a recurrence state. 

Proof: Consider a system configuration a with u min < l max - Upon addition of 2 particles to sites 1 
and 2N—1 and 1 particle to sites 2 and 2N, it is easy to verify that exactly one toppling will occur 
in the following sites 1,2,..., 2u min -2, 2u min , 2l max -l, 2l max +l, 2l max +2, . . . , 2N. After that, the 
avalanche stops because both sites 2u min — 1 and 2l max "absorb" the incoming particles and prevent 
further toppling. Since u min < l max , sites such as 2u min — 1 and 2l max will not topple during the 
avalanche and hence after the system will not relax back to a. Thus a (jL fi. □ 

Thus the remaining possible recurrence state configurations are those a = (hi) with —3 < 
hn-\ < and — 1 < h%i < for i = 1, 2, . . . , N, and with u m i n > l ma x- The total number of such 
states, T, is given by: 



TV 



A? 



EE 2 ' 

i=l j=l 



-1 2 A N~ 



l 2 r 



E 2i 



-i 2 4 iv- 



+ 2 N 2 N 



(10) 



where the first term is the number of configurations with u m i n < N and l max > 1, second term is 
the number of configurations with l max = 0, and third term is the number of configurations with 
Umin — N+l. After some computation, we find T = (N + 1)4 N = det A, and hence the set of all 
recurrence states of the system is 



n 



{a = 
for % 



: a 2i _! G {-3, -2, -1, 0}, a 2j E {-1, 0} 



1,2, 



,N and u min (a) > l min 



»}■ 



(11) 



Now, we go on to show that the system is indeed self-organized critical. Unlike Example 1, it is 
not easy to argue the distribution of avalanche size D(s) owing to the complexity of the recurrence 
phase space Q. So we take the alternative approach by calculating the two point correlation function 
Gij of the system which is defined as the average number of toppling occurs in site j given that a 
particle is introduced to site i, and is given by = A^ 1 0. 

In the Appendix, we show that 



Co,; 



2i,2j 



j{N-i) 
4(AT+1) 

(i+l)(N+l-j) 
4(JV+1) 



if % > j 
otherwise 



12a) 



G 



2i,2j-l 



4(7V+1) 

(i+l)(N+l-j) 
4(7V+1) 



if i + 1 > j 
otherwise 



(12b) 



G 



2t-l,2j 



(i-l)(N+l-j) 
2(N+1) 

j(N+2-i) 
2(N+1) 



if i - 1 < j 
otherwise 



(12c) 
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and 



G 2i . 12j -i={ 2(N+1) , (12d) 

^~2(jv+i)~^ otherwise 

for i,j = 1,2, ... ,N. Obviously, the two point correlation function varies linearly with the 
distance between sites i and j. Upon a uniform and random particle addition, 1/f 2 scaling is 
observed as N—>oo and the model is indeed self-organized critical (but with a trivial exponent). 

This model is not equivalent to any "ordinary" one-dimensional ASM with toppling rules being 
translational invariant except possibly at the boundary. Suppose the contrary, we can find a set of 
h™ ew = fi(hi, . . . , h 2 N) such that the system becomes a one-dimensional ASM after applying these 
transformations. f\, however, is independent of hj for all j '• ^ i or else toppling of site i must depend 
on the neighboring sites, making the transformed system not an ASM. It is easy to check that any 
transformation fi = fi(hi) cannot reduce Eq. (^) to a Toeplitz form and hence this model is not 
equivalent to any "ordinary" one- dimensional ASM. 



IV. CONCLUSIONS 

In summary, we have explicitly constructed the simplest possible class of cellular automaton 
examples exhibiting self-organized criticality without the presence of local particle conservation law: 
namely, Abelian Sandpile Models on one-dimensional superlattices. Moreover, a simple method 
of finding recurrence phase space configurations, based on the idea of the burning algorithm is 
introduced, which is useful when the dimension of the system is low. 



APPENDIX: APPENDIX: FINDING A 

We rearrange the site labels of the system using the map: 

2% — > i 



-i 



2i-l — N + i (A1) 



for i = 1, 2, . . . , N, the toppling matrix A can be re- written as 

(A2) 



A = A {2N) 

"new — ^new 



2D T B 
2B T D 



where B T denotes the transpose of B. Moreover, B and D are Ax A-matrices whose elements are 
given by 

^ I otherwise (A3a) 



and 



2 ifi = j 

Dij = { -1 if i-3 = 1 (A3b) 
otherwise 
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respectively. According to the block matrix inversion formula 

A 



-i 

'new 



0.5 X- 1 

-D- 1 B T X- 1 



0.5 {D T Y'BY- 1 
Y- 1 



(A4) 



where X = D T — B D l B T and Y = D — B T (D T ) 1 B. Since D is a Toeplitz matrix, D 1 can be 



evaluated easily [13 and is given by 



D 



-i 



0.5^' +1 




if i > j 
otherwise 



(A5) 



As a result, 



X, 



( 2 

-1 
1.5 

-0.5^' +1 




if i = j = N 

if i ~ j = 1 

if % = j and i < N 

if i > j and i < N 

otherwise 



(A6) 



We divide X into four block matrices by partitioning the sites into two sets, namely: {1, 2, ... , N—l} 
and {iV}. The only (N — l) x (N — l)-matrix so formed is in Toeplitz form whose inverse can be 
found readily |13| . After that, by means of the block matrix inversion formula again, we obtain 



X 



j(N-i) 
2(7V+1) 

(i+l)(AT+l-j) 
2(7V+1) 



Using the same method, we find that 



Y 



-i 



j(N+2-i) 
2{N+1) 

(i-l)(N+l-j) 
2(N+l) 



if i > j 
otherwise 

if % > j 
otherwise 



(AT) 



(A8) 



Now D- 1 B T X- 1 and (D T ) 1 B Y~ x can be evaluated using Eqs. (^|), (fA^) and flggp . Thus all 
four block matrices in Eq. (\K/§) are computed. Finally, Eqs. (^2]a)-([l2|d) are obtained by changing 
the labels of the sites back to the original ones. 
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